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Abstract 

BEAM  (Beam  Element  Analysis  with  Mechanisms)  is  a  finite  element  program  for  the 
analysts  of  the  collapse  of  vehicle  structures.  BEAM  can  quickly  estimate  the  collapse  load 
and  plastic  hinge  locations  in  three-dimensional  thtn-walled  frames  that  would  collapse  in 
bending.  Such  information  can  be  combined  with  separate  rigid-plastic  modelling  to  enable 
the  large-deflection  design  of  vehicle  frameworks  so  that  their  collapse  in  roll-over 
accidents,  for  instance,  can  be  prevented.  One  advantage  of  BEAM  is  that  it  allows  a  user 
to  define  his  own  experimentally  or  analytically  determined  loads  and  moments  as  failure 
criteria.  The  theory  and  operation  of  BEAM  are  presented  and  its  application  to  the 
modelling  of  a  collapsing  two-dimensional  frame  is  described. 
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BEAM:  A  Finite  Element  Program  for  the 
Collapse  Analysis  of  Vehicle  Structures 


1.  Introduction 

Accurate  large-deflection  design  techniques  have  been  developed  for  vehicles 
whose  superstructures  are  constructed  of  thin-walled  frames  which  fail  by 
collapsing  in  bending  and  exhibit  drooping  load-deflection  behaviour  (1),(2).  The 
impetus  for  the  development  of  these  methods  has  come  from  the  requirement  to 
design  passenger  buses  to  meet  new  Australian  Design  Rules  for  occupant 
protection  in  roll-over  accidents  (3).  These  same  techniques  can  also  be  applied  to 
military  vehicles  which  use  their  ippier  body-work  as  a  Roll-Over  Protective 
Structure  (ROPS)  or  even  for  the  design  of  add-on  roll-bars. 

The  entire  load -deflection  behaviour  of  thin-walled  vehicle  framewoi  ks  can  be 
determined  by  simply  superimposing  separate  elastic  and  rigid-plastic  models  of 
a  frame's  behaviour  (l)/(2).  While  in  many  cases  vehicle  frames  can  be  designed  to 
withstand  roll-over  accidents  with  two-dimensional  techniques,  there  are  some 
instances  where  it  is  important  to  include  the  effects  of  structural  members  that 
are  orientated  in  a  third  dimerrsion. 

A  finite  element  technique  was  considered  to  be  the  most  sensible  method  for 
determining  the  influence  of  three-dimensional  effects  on  frame  collapse 
behaviour.  The  three-dimensional  finite  element  program  'BEAM'  (Beam  Element 
Analysis  with  Mechanisms)  (2)  was  specifically  written  for  this  purpose  because 
conventional  finite  element  programs  were  unsuitable.  One  advantage  of  BEAM 
over  other  programs  is  that  it  allows  a  user  to  define  their  own  experimentally  or 
analytically  determined  loads  and  moments  as  failure  criteria.  BEAM  will  allow  a 
quick  determination  of  the  collapse  load  of  a  vehicle  frame  as  well  as  its 
deformation  up  to  this  load,  in  the  process  predicting  the  locations  of  the  plastic 
hinges  which  form  its  collapse  mechanism.  The  latter  mformation  allows  a  rigid- 
plastic  modelling  technique  to  be  easily  applied  to  determine  the  behaviour  of  the 
vehicle  frame  beyond  the  point  of  collapse. 

This  report  outlines  the  theory  and  operation  of  BEAM  and  its  application  to 
modelling  a  collapsing  thin-walled  frame. 


2.  Finite  Element  Modelling  Collapsing  Thin- 
Walled  Vehicle  Structures 

For  various  reasons  it  can  be  difficult  to  model  collapsing  vehicle  frames  with  the 
finite  element  method.  This  is  partly  because  of  the  gross  plastic  deformation  and 
geometric  distortion  associated  with  the  formation  of  plastic  hinges  within  the 
thin-walled  sections  of  the  vehicle  frame.  Specialised  large-deformation  finite 
element  programs  are  used  to  model  this  deformation  because  they  have  specific 
algorithms  for  accurate  and  time-efficient  modelling  of  gross  plastic  deformation 
and  distortion,  e.g.  (5).  General  purpose  nonlinear  finite  element  programs,  that 
are  used  for  common  engineering  structural  analyses,  for  instance,  are  too 
inefficient  for  these  purposes  (6),  (7). 

Specialised  large-deformation  programs  have  major  limitations.  These  programs 
may  require  access  to  expensive  computer  systems  and  require  considerable 
specialised  manpower  and  time  to  code  and  run  a  vehicle  collapse  analysis.  In 
addition  to  these  efficiency  aspects  there  are  also  important  technical  problems 
which  have  not  yet  been  overcome. 

Even  the  largest  and  most  specialised  large-deformation  programs  have  not  yet 
been  proven  to  accurately  predict  fundamental  failure  mechanisms  (8).  While  such 
programs  can  still  produce  excellent  representations  of  experimentally  deformed 
components,  quantitative  comparisons  of  numerical  and  experimental  crush  force 
are  generally  not  available  (8).  In  any  case,  the  accuracy  of  such  comparisons,  if 
available,  can  be  greatly  changed  by  the  values  assumed  for  various  ill-defined 
quantities,  for  instance,  the  strain  rate  effect  on  flow  stress  (8).  In  some  cases  the 
manual  intervention  to  the  modelling  process  that  is  required  to  produce  accurate 
failure  mechanisms  can  have  little  relationship  to  reality',  e  g.  (V).  To  obtai,'. 
realistic  results  from  such  numerical  simulations,  it  appears  that  the  deformed 
geometry  of  a  structure  is  required  a  priori.  This  would  require  a  vehicle  prototype 
to  be  built  and  crash  tested,  an  unsatisfactory  situation. 

Furthermore,  the  large  computer  systems  required  to  efficiently  run  large- 
deformation  programs  may  not  provide  the  continuous  access  needed  for  design 
improvement  and  optimisation  in  the  early  stages  of  product  development  when 
the  base  structural  design  is  being  finalised.  Precise  modelling  is  unnecessary,  as 
at  this  stage  the  exact  design  details  are  generally  not  specified.  What  is  needed  at 
t.he  initial  design  stages  is  a  quick  first-order  approximation  of  the  structural 
strength  of  the  major  components.  Such  analyses  can  be  used  to  optimise  the  base 
structural  strength  of  a  design  well  before  vehicle  prototypes  have  to  be  built.  A 
vehicle  model  can  therefore  be  simplified  to  just  the  major  load  bearing  members, 
for  instance,  a  framework  consisting  of  beam  elements. 

A  promising  improvement  to  modelling  gross  plastic  deformation,  appears  to 
be  to  incorporate  heuristic  or  hybrid  elements  into  a  conventional  finite  element 
mesh  (8).  This  approach  has  been  adopted  and  developed  by  workers  such  as 
Kecman  (10),(11),  Mahmood  and  Paluszny  (12)  and  Wierzbicki  and  Abramovich 
(13).  Heuristic  type  elements  can  incorporate  the  mechanics  of  the  plastic  hinges 
found  in  the  grossly  deformed  regions  of  thin-walled  sections  (8).  The  use  of  such 
elements  can  eliminate  many  thousands  of  conventional  elements  from  an  analysis 
as  the  plastic  hinges  are  simply  treated  as  'black  boxes'.  In  other  words,  the  exact 
details  of  the  deformation  of  the  plastic  hinges  are  not  important  to  the  operation 
of  the  program,  only  their  external  effects,  i.e.  their  externally  recognised  resisting 
moments.  The  resulting  coarser  numerical  models  require  less  hardware  and 
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human  resources,  and  allow  an  analysis  of  a  crash  to  be  performed  on  a  PC  or 
workstation. 

There  are  various  ways  of  incorporating  heuristic  elements  into  finite  element 
programs.  The  'CRASH-D'  program  (CRAnfield  Structural  High  Deformation)  or 
as  it  is  now  known,  'CIC-STAT',  is  the  best  known  example  of  this  approach  (14). 
This  program  uses  beam  elements  to  analyse  the  structural  collapse  of  three- 
dimensional  frames  to  large  deflections.  The  resistance  of  the  plastic  hinges  in  a 
collapsing  frame  is  simply  obtained  from  analytical  or  expierimentally  determined 
quasi-static  moment-rotation  curves,  allowing  the  specific  details  of  the 
deformation  at  the  hinges  to  be  ignored.  Consequently  a  crash  analysis  can  be 
much  more  quickly  pci  formed  than  if  conventional  finite  element  modelling 
methods  were  used.  Applications  have  included  the  crushing  of  various  vehicles 
and  components,  viz,  truck  cabins  (15),  door  intrusions  and  roof  crushes  of  cars 
(16),(17),(18)  as  well  as  bus  roll-overs  (19),(20). 

Programs  like  CRASH-D  can  be  used  to  model  the  collapse  of  vehicle  structures 
in  crashes  that  are  amenable  to  static  analyses,  with  dynamic  effects  being  ignored 
(21 )  Passenger  bus  roll-over  accidents,  for  instance,  fall  into  this  category.  This  is 
in  part  because  such  accidents  tend  to  occur  laterally  at  low  speeds  (20), (22). 
Inertial  effects  can  effectively  be  neglected  because  deformation  is  largely 
confined  to  the  superstructure  of  the  buses,  with  most  of  their  mass  being 
concentrated  at  floor  level.  In  addition,  it  is  often  the  case  that  the  dynamic 
collapse  modes  in  bending  are  the  same  as  the  static  ones  (11).  At  least  this  is  the 
case  for  the  square  thin-walled  tubes  from  which  passenger  buses  have  been 
constructed.  More  importantly,  however,  the  dynamic  enhancement  of  the 
material  yield  stress  as  well  as  flow  stress  at  large  plastic  strains  can  be  accounted 
for  to  some  extent  by  multiplicative  factors  (23),(24). 

Program  BEAM  can  be  applied  to  the  same  class  of  problems  as  CRASH-D  and 
accounts  for  a  plastic  hinge  in  a  member  by  modifying  the  element  stiffness  matrix 
to  include  a  pin-joint  which  will  represent  the  physical  effect  of  a  plastic  hinge,  i.e. 
zero  stiffness.  This  operation  is  performed  when  an  experimentally  determined 
failure  criteria  is  exceeded,  the  pin-joint  being  capable  of  modelling  the  behaviour 
of  a  plastic  hinge  from  its  point  of  formation  up  to  when  drooping  moment- 
rotation  behaviour  is  exhibited.  Therefore  BEAM,  unlike  CRASH-D,  cannot 
determine  the  drooping  load-deflection  curve  of  a  collapsing  thin-walled  frame. 
However,  the  load-deflection  curve  that  BEAM  calculates  can  be  combined  with  a 
separate  rigid-plastic  analysis  of  the  collapsing  frame.  This  allows  the  behaviour 
of  such  a  frame  to  be  dete'mined  to  large  deflections  where  drooping  load- 
deflection  behaviour  is  exhibited. 


3.  Program  BEAM 

3.1  Step-By-Step  Analysis  and  Program  BEAM 

Simple  Plastic  Analysis  has  been  shown  to  be  capable  of  accurately  predicting  the 
collapse  load  of  simple  thin-walled  rectangular  frames  (25).  This  method  requires 
an  accurate  assumption  of  the  collapse  mechanism  of  a  structure,  which  is 
straightforward  for  simple  structures.  However,  for  more  complex  structures. 
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trial  and  error  procedures  involving  significant  intuition  are  required  to 
determine  the  correct  collapse  mechanism.  It  is  difficult  to  program  computers  to 
make  the  necessary  judgements  for  these  procedures,  but  computer  methods 
using  Step-By-Step  analysis,  for  instance,  have  been  develop>ed  that  can  quickly 
calculate  the  collapse  load  and  plastic  hinge  locations  of  statically  loaded  frames. 
Wemg  (26)  first  used  this  method  to  determine  the  collapse  load  of  two- 
dimensional  frames. 

Step-By-Step  analysis  involves  making  an  initial  elastic  analysis  of  a  frame 
under  a  working  load  and  then  dividing  the  moments  at  the  end  of  each  element 
into  their  full-plastic  values  to  obtain  load-factors  to  form  plastic  hinges;  the 
smallest  load-factor  identifying  where  the  first  plastic  hinge  would  form.  The 
relevant  element  stiffness  matrix  is  then  modified  to  include  a  pin-joint  at  the 
appropriate  node  by  zeroing  the  relevant  stiffnesses.  The  process  is  repeated  and 
successive  elastic  analyses  performed  until  the  structure  is  converted  into  a 
mechanism. 

Program  BEAM  was  written  in  FORTRAN  77  as  a  three-dimensional  beam 
element  finite  element  program.  It  is  modular  in  nature  and  employs  an  upper- 
half-bandwidth  solver  that  uses  Gaussian  Elimination  and  Back  Substitution  to 
obtain  a  solution  (27).  The  Step-By-Step  procedure  of  BEAM  was  originally 
derived  from  the  two-dimensional  EPFOSS  program  of  Harrison  (28).  This  latter 
program  originally  accounted  for  bending  and  axial  deformation  only  and  used  a 
Triple  Decomposition  method  for  solution.  The  12  x  12  linear  elastic  stiffness 
matrix  of  BEAM  was  able  to  account  for  bending,  shear,  axial  and  torsional 
deformation.  Fig.  1  defines  the  degrees  of  freedom  of  this  element  However,  the 
current  version  of  BEAM  only  checks  for  the  possibility  of  bending  collapse  in  the 
x  and  y  local  axes  of  a  section.  Furthermore,  BEAM  currently  has  no  provision  to 
detect  plastic  hinges  formed  under  significant  multi-axial  loads.  This  limited 
capability  is  sufficient  for  the  collapse  of  regular  prismatic  frames,  for  instance 
civilian  passenger  buses,  but  other  failure  modes  can  contribute  to  the  structural 
collapse  of  frames  of  more  general  geometry.  Axial  crushing,  e  g  (29)  and 
torsional  collapse  should  also  be  accounted  for.  Further  experimental 
characterisation  is  also  required  to  determine  correctly  the  failure  criteria  for  these 
modes  of  collapse  as  well  as  in  multi-axial  situations.  BEAM  has  provision  to 
include  these  failure  criteria  when  appropriate  data  becomes  available,  e  g.  (30) 

BEAM  employs  a  linear  elastic  stiffness  matrix  that  is  based  on  the  original 
undeformed  structure  geometry.  This  is  because  the  aim  of  BEAM  is  to  provide 
only  a  quick  first-order  estimate  of  a  frame's  load-displacement  curve  up  until  it 
forms  a  mechanism.  Consequently  BEAM  cannot  account  for  effects  such  as  the 
P-8  Effect  (4),  i  e  when  a  strut  carrying  an  axial  load  P  deflects  a  lateral  distance  8, 
its  bending  stresses  are  increased.  Nor  can  BEAM  account  tor  the  reduction  of 
plastic  moment  capacity  due  to  axial  loads  (4)  While  non-lmear  elastic  stiffness 
matrices  can  be  included  that  can  account  for  a  structure  s  changing  geometry, 
they  would  make  the  operation  of  BEAM  cumbersome  and  detract  from  its 
inherent  simplicity.  The  load-deflection  response  of  a  real  structure  w  ill  therefore 
not  be  as  stiff  as  that  predicted  by  BEAM.  However,  first-order  analysis  is 
sufficiently  accurate  for  quick  analyses  of  thin-walled  structures  up  tc  -rollapse 
The  errors  involved  with  such  a  philosophy  are  likely  to  be  minimal  considering 
the  large  deformations  that  such  structures  experience.  Furthermore,  if  is  assumed 
that  a  structure's  joints  and  connections  maintain  their  integrity  to  large 
deformations.  Since  it  is  not  easy  to  predict  joint  strength  to  large  deflections,  this 
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should  ideally  be  verified  experimenfallv  prior  to  nruxlelling  by  either  individual 
joint  tests  (31)  or  full-sized  frame  tests  (2) 


Figure  1  Degrees  Freedom  ussoemted  unth  the  three-dimensional  beam  element  of 
BF.AM  that  is  defined  f’v  nodes  i  arid  ;  Arrows  with  single  heads  denote  loads  while 
double-headed  arrows  denote  moments 


The  insertion  of  a  plastie  hinge  (or  pin-joinf)  can  be  achieved  bv  modifying 
portions  of  the  relevant  element  stiftncss  matrm.  As  an  example.  Fig  2  shows 
possible  element  stiffness  matrix  modifications  for  a  simple  beam  element  with 
two  bending  degrees  of  freedom  The  bending  stiffness  at  the  nixie  of  the  plastic 
hinge  and  the  cross-stiffnesses  were  zeroed,  witn  the  bending  stiffness  at  the 
opposite  end  of  the  element  being  reduced  to  three-quarters  of  its  usual  stiffness 
As  has  been  the  practice  with  other  workers  (14),(26),(28),  the  axial  stiffnesses 
were  not  modified.  However,  the  12  x  12  stiffness  matrix  of  BEAM  also  includes 
shear  stiffness  terms  which  must  be  modified  to  account  for  the  presence  of 
plastic  hinges.  The  appropriate  shear  stiffnesses  were  modified  according  to 
established  rule?  (32),  However,  it  was  assumed  {for  simplicity)  that  a  plastic 
hinge  forming  on  one  kxal  element  axis  would  not  influence  the  bending 
stiffnesses  associated  with  the  other  local  element  axis  Moreover,  the  formation 
of  a  plastic  hinge  in  bending  was  not  considered  to  affect  a  member  s  torsional 
stiffness.  While  these  assumptions  may  not  be  strictly  correct,  they  are 
nevertheless  reasonable  in  the  context  of  the  present  sparse  knowledge  of  the 
three-dimensional  properties  of  plastic  hinges  in  bending  collapse. 


Plastic  hinge  at 
end  i 


0 - O 

'  Plastic  hinge  at  both  i 
ends 


Figure  2  Element  stiffness  matrix  modifications  caused  fry  iiarwus  plastic  hinge 
insertions  for  a  simple  beam  element  with  two  bending  degrees  of  freedom  (After  Harrison 


(28)). 


There  are  other  assumptions  inherent  in  the  Step-By-Step  priKcdure  of  BEAM 
that  are  common  to  Simple  Plastic  Analysts  (23)  For  instance,  it  is  assumed  that  as 
the  load-factor  is  increased  until  the  structure  reaches  its  ultimate  load,  each 
plastic  hinge  maintains  a  constant  moment  It  can  also  sometimes  ivcur  that  a 
plastic  hinge  that  forms  early  in  an  analysis  may  not  be  required  for  the  final 
frame  collapse  mechanism  The  moment  at  such  a  hinge  can  actually  unload 
While  this  phenomena  cannot  be  accounted  for  in  BEAM  (28),  the  predicted  load 
factor  for  collapse  would  nevertheless  be  conservative,  i  e  a  lower-bound 
estimate,  because  the  equilibrium  and  yield  conditions  are  satisfied  but  not  the 
mechanism  condition  (4)  This  is  an  example  of  a  partial  collapse  (33) 

One  special  feature  of  BEAM  is  the  FORMRT  subroutine  that  calculates  the 
three-dimensional  rotation  matrices  that  are  used  to  transform  element  stiffness 
matrices  into  global  co-ordinates  FOR.MRT  was  based  on  the  vector  method  of 
Wilson  (34),  and  allow's  BEAM  to  be  easily  applied  to  three-d.mensional  frames 
where  the  orientation  of  the  principal  axis  (i.e.  plane  of  principal  bending)  has  to 
be  specified.  Three  nodes  are  required  to  be  specified,  viz,  nodes  i  and  )  defining 
an  element's  location  and  nixie  k,  that  defines  its  principal  axis.  Fig  3  The 
program  requires  nixie  k  to  be  specified  in  the  positive  quadrant  of  the  local  axes 
defined  by  the  first  and  second  degrees  of  freedom  of  the  element 


X 

Figure  .1  Node  k  defines  the  principal  axis  of  the  element  defined  by  nodes  i  and  j.  Node  k. 
denoted  by  (X),  must  be  specified  within  the  positive  quadrant  defined  by  the  axes  of  the  first  and 
second  degrees  of  freedom  of  the  element  (After  Wilson  (34)). 

3.2  Operation  of  Program  BEAM 

The  operation  of  program  BEAM  is  summarised  in  the  flow  chart  in  Fig.  4.  A  full  listing 
of  the  source  code  for  BEAM  is  given  in  a  separate  publication  (35).  Appendix  1  outlines 
the  layout  of  a  generic  input  file.  BEAM  gives  an  analyst  the  choice  of  either  an  Elastic 
analysis  or  a  Step-By-Step  analysis  to  collapse.  If  the  latter  option  is  chosen,  the 
structural  stiffness  matrix  (GSTIF)  is  assembled  in  subroutine  ASSEMB  which  obtains 
the  elemental  stiffness  matrices  from  subroutine  EBEAM  and  the  corresp>onding  rotation 
nruitrices  from  subroutine  FORMRT.  Matrix  (GSTIF)  is  then  inspected  for  zero  stiffnesses 
in  its  main  diagonal,  a  singular  matrix  indicating  that  the  structure  has  been  transformed 
into  a  mechanism.  This  check  is  the  pnmary  means  for  terminating  the  program.  If  the 
analysis  continues,  a  typical  working  load  is  assigned  to  the  loading  matrix  ( AJJ.  Matrix 
(CiSTlF)  is  then  corrected  for  restrained  boundary  conditions  in  subroutine  BOUND  and 
solved  by  subroutine  SOLVE  to  obtain  the  nodal  displacements. 

BEAM  incorporates  another  test  to  check  for  collapse  because  round-off  errors  could 
delay  program  termination  if  a  check  was  made  only  for  singularity.  An  additional 
check  is  made  for  very  large  computed  displacements.  l,arge  displacements  indicate  that 
the  structure’s  load-displacement  curve  had  effectively  become  horizontal,  as  would  be 
the  case  in  the  early  stages  of  structural  collapse. 


13 


The  member  stress  resultants  matrix  [WRKSRl  and  member  displacements  matrix 
[DM]  are  next  calculated.  The  moments  at  the  bending  degrees  of  freedom  are  then 
divided  into  the  specified  failure  moments  to  obtain  load  factors  for  possible  plastic 
hinges.  The  first  plastic  hinge  will  form  at  the  degree  of  freedom  with  the  smallest  load 
factor,  SLF,  BEAM  determining  whether  this  hinge  forms  in  the  local  member  y  axis  or  z 
axis.  This  information  is  recorded  along  with  which  end  of  the  element  the  plastic  hinge 
is  situated.  Matrices  [WRKSR],  [DM]  and  the  nodal  displacements  are  therefore 
multiplied  by  SLF  and  added  to  the  cumulative  stress-resultants,  matrix  [CUMSRj; 
cumulative  member  displacements,  matrix  [CUMDM];  and  the  cumulative  nodal 
displacements,  matrix  [CUMX];  respectively.  SLF  is  also  added  to  CUMLF,  the 
cumulative  load-factor.  These  quantities  are  printed  together  with  the  location  and  load- 
factor  of  the  plastic  hinge. 

The  program  loop  is  repeated  with  the  structure  stiffness  matrix  [GSTIFJ  again  being 
assembled.  Once  plastic  hinges  have  begun  to  form,  the  element  stiffness  matrices  that 
are  generated  within  subroutine  EBEAM  are  modified  to  allow  for  the  presence  of 
plastic  hinges  in  either  or  both  ends  of  the  element  on  its  y  or  z  local  member  axes.  These 
heuristically  modified  matrices  are  again  assembled  in  ASSEMB  and  the  kxip  continued 
with  the  original  loading  matrix  ( AJl  again  applied  and  a  further  solution  i  btained. 

Eventually,  when  sufficient  plastic  hinges  have  formed  to  initiate  structural  collapse, 
the  program  will  detect  either  a  zero  diagonal  in  the  structure  stiffness  matrix  or 
excessive  nodal  displacements.  The  program  then  proceeds  to  line  490  and  sets  a  flag  so 
that  a  final  output  Cim  be  made  of  matrices  [CUMX|,  (CL'.MSK)  and  [CUM DM)  along 
with  the  unfactored  displacements  remaining  after  the  last  solution. 
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An  important  feature  of  BEAM  is  its  ability  to  insert  plastic  hinges  of  zero 
stiffness  into  the  structure  stiffness  matrix  without  creating  a  singular  matrix.  This 
is  achieved  because  each  node  has  some  finite  stiffness  associated  with  it,  that  is  if 
a  plastic  hinge  forms  at  one  end  of  an  element,  the  adjacent  element  connecting  to 
that  element  will  still  have  some  finite  stiffness.  This  allows  the  lack  of  stiffness 
associated  with  a  plastic  hinge  to  be  incorporated  into  an  elastic  analysis.  One 
problem  that  was  encountered  during  validation  was  the  inability  of  BEAM  to 
insert  a  plastic  hinge  into  a  degree  of  freedom  associated  with  a  fixed  boundary 
condition  that  had  zero  displacement.  BEAM  overcame  this  problem  by 
incorporating  a  test  in  subroutine  EBEAM  which  checked  whether  any  plastic 
hinges  coincided  with  such  fixed  boundary  conditions  If  coincidence  was 
detected,  the  appropriate  rotation  degree  of  freedom  was  set  to  1  x  10'^^,  a 
negligible  value  but  sufficient  to  prevent  the  structure  stiffness  matrix  from 
becoming  singular 

In  comparison  to  CRASH-D  (14)  BE.AM  does  not  appear  to  account  for  the 
resisting  moments  of  existing  plastic  hinges  in  its  solution  procedure,  the  solution 
technique  of  BEAM  being  purely  stiffness  based.  However,  the  load-displacement 
curve  of  BEAM  is  obtained  from  the  cumulative  stress-resultants  matrix 
[CUMSR].  The  moments  in  this  matrix  at  the  degrees  of  freedom  associated  with 
plastic  hinges  are  equal  to  the  failure  moment  because  the  results  of  each  elastic 
analysis,  matrix  [WRKSR],  are  multiplied  by  the  minimum  load  factor,  SEE,  before 
being  added  to  matrix  [CLMSR]  (28), (.4h)  Cine  advantage  of  BE.AM  is  its  abilitv  to 
quickly  determine  a  structure  s  load-displacement  curve  up  to  the  collapse  load 
with  a  minimum  number  of  solution  steps  In  comparison,  the  incremental 
solution  procedure  of  CRASH-D  requires  more  operator  skill  and  an  indefinite 
number  of  structure  stiffness  matrix  inversions. 


3.3  Verification  of  Program  BEAM 

BEAM  was  elastically  verified  in  three  dimensions  by  analysing  the  frame  shown 
in  Fig.  5. 

The  Step-By-Step  procedure  of  BEAM  was  verified  by  analysing  the  frames 
detailed  in  Fig.  6  and  Fig.  7  up  to  the  point  of  collapse.  These  frames  had  1  and  6 
degrees  of  redundancy,  respectively,  and  so  required  only  2  and  7  plastic  hinges, 
respectively,  to  form  for  structural  collapse.  As  well  as  predicting  the  correct 
locations  and  order  of  plastic  hinge  formation,  BEAM  obtained  exactly  the  sanr.e 
displacements  and  loads  as  the  two-dimensional  EPFOSS  program  of  Harrison 
(28).  The  two-dimensional  SOAPS  program  (37),  which  used  a  similar  Step-By- 
Step  method,  independently  confirmed  these  results.  The  two-dimensional  PCAP 
program  (38),  which  used  a  linear  programming  technique  for  solution,  again 
independently  confirmed  the  collapse  load  of  these  frames.  Overall,  the  operation 
of  BEAM  was  sufficiently  verified  to  justify  its  application  to  predict  the  plastic 
hinge  locations  and  collapse  load  of  an  experimentally  tested  two-dimensional 
frame. 


4) 
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(a) 

Translations  (mm)  Rotations  (deg) 

X  Y  Z  X  Y  Z 

1  0  0  0  0  0  0 

2  0  0  0  0  0  0 

3  0  0  0.147  -0.0267  0.1274  0 

4  0  0  9.117  -0.0267  0  0 


(b) 


Figure  5;  Frame  1,  (a)  initial  loading  and  (b)  results  (global  nodal  displacements).  ^ 

Member  3  is  310UC97  with  a  vertical  web.  Members  1  and  2  are  3WUB46  with  the  webs 

lying  in  the  plane  containing  the  relevant  member  itself  and  Member  3.  Young's 

Modulus,  E,  and  the  Shear  Modulus,  C,  are  206x10^  kNm  ^,  and  90x10^  kNm'^, 

respectively. 

• 


19 


•  • 


(b) 

Figure  7:  Frame  3:  (a)  initial  geometry,  showing  node  numbers  with  element 
designations  circled;  and  (b)  deformed  geometry  with  sequence  of  plastic  hinge  formation. 
(After  Harrison  (28)). 


4.  Modelling  the  Collapse  of  a  Two- 
Dimensional  Frame  with  Program  BEAM 

A  two  dimensional  frame  made  from  50x50  mm  square  thin-walled  tubes  of  2  mm 
wall  thickness.  Frame  4,  was  tested  to  confirm  the  validity  of  the  model 
predictions  (2).  Fig.  8  depicts  the  geometry  of  the  frame,  (a),  and  its  idealised  finite 
element  model,  (b).  The  experimental  load-displacement  response  of  the  frame 
was  obtained  by  simply  displacing  the  frame  by  means  of  a  hydraulic  tension  jack 
that  pulled  from  an  initial  angle  of  15°  (2).  The  applied  load  and  tlie  movement  of 
the  point  of  load  application  were  measured  as  deformation  proceeded. 
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Figure  8:  The  initial  geometry,  (a),  and  finite  element  model,  (b).  of  the  tested  hoo- 
dimensional  thin-walled  frame,  Frame  4.  Structural  nodes  in  <b)  are  depicted  by  {•),  k 
nodes  (refer  to  Fig.  3)  by  (X),  and  the  element  numbers  circled. 


•  • 


An  input  file  was  prepared  with  a  working  load  of  5  kN,  The  input  file  and  the 
output  file  are  included  in  Appendix  2  and  Appendix  3,  respectively.  Pure  ^ 

bending  experiments  were  used  to  determine  the  maximum  moment  of  the  tubing 
(2),  (39).  A  maximum  moment  of  3.12  kNm  was  obtained  and  u.sed  for  the  failure 
moment  in  the  input  file.  Note  that  in  the  output  file  negative  signs  denote  the 
locations  of  the  plastic  hinges.  The  plastic  hinges  formed  first  at  the  base  and  then, 
with  further  deformation,  at  the  roof  joints  as  expected.  The  plastic  hinges  formed 

in  the  following  order;  • 

plastic  hinge  1  at  node  1  at  a  load  of  10.466  kN 
plastic  hinge  2  at  node  4  at  a  load  of  10.484  kN 

plastic  hinge  3  at  node  2  at  a  load  of  12.173  kN  * 

plastic  hinge  4  at  node  3  at  a  load  of  12.176  kN 
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Horizontal  load,  kN 


These  Step-By-Step  analysis  results  were  confirmed  independently  usinj;  the  Step-By- 
Step  SOAPS  (37)  and  PCAP  (38)  finite  element  programs  and  are  plotted  as  a  function  of 
the  horizontal  displacement  of  the  loading  point  in  I'ig.  9.  Separate  clastic  and  rigid- 
plastic  mixlels  are  also  superimposed.  It  has  already  been  shown  that  the  complete  load- 
displacement  curve  of  this  frame  can  be  mcxlelled  accurately  by  superimposing  separate 
elastic  and  rigid-plastic  models  of  its  behaviour  (1),(2).  This  was  for  a  first-order  elastic 
model  which  assumed  that  the  frame  behaved  elastically  up  to  its  collapse  load  and  fi>r  a 
rigid-plastic  model  which  assumed  that  the  plastic  hinges  nece.ssary  for  the  collapse  of 
the  frame  began  rotating  at  zero  displacement. 


Horizontal  displacement,  mm 


Figure  9:  The  load-displacement  curve  for  Frame  4,  comparing  elastic,  rigid  plastic  and  titqi 
Bp-Step  analyses  with  experimental  results.  The  combined  model  of  the  load  displacement  curve 
is  shmvn  in  hold.  The  rigid-plastic  model  assumes  that  the  plastic  hinges  form  at  the  corners  of 
the  frame  and  its  boundaries  with  the  fixed  supports. 


I'igure  9  shows  that  if  the  results  of  the  Step-By-Step  analysis  are  combined  with  the 
separate  rigid-plastic  analysis  then  a  better  model  of  the  complete  load-displacement 
curve  of  the  frame  can  be  obtained,  as  shown  by  the  bold  lines.  Fig.  10  has  an  expanded 
displacement  scale  to  show  more  clearly  the  detail  of  the  first  120  mm  displacement.  A 
line  between  the  Step-By-Step  analysis  and  the  collapsing  portion  of  the  rigid-plashc 
curve  can  be  sketched  by  hand  to  complete  the  mcKlel  (1),(2)  (shown  by  the  bold  dashed 
lines  in  Figures  9  and  10). 

It  was  fi>und  in  the  experimental  frame  test  that  a  maximum  load  of  12..30  kN 
was  reacht“d  at  88  mm  deflection,  with  collapse  cKCurring  at  a  load  of  12.20  kN  at 
101  mm  (.leflection.  BF.AM  actually  determines  the  point  at  which  the  frame  first 
forms  a  mechanism.  At  this  point  BEAM  prtxlicted  a  maximum  loatl  of  12.18  kN 
for  ,39  mm  deflection.  Once  plastic  hinges  form  in  the  particular  thin-walleil 


members  of  the  frame,  they  deform  at  a  relatively  constarrt  moment  before 
exhibiting  a  drcxjping  moment-rotation  curve  (39).  As  a  consequence,  a  collapsing 
thin-walled  frame  will  also  deform  at  a  relatively  constant  load  from  the  point  of 
mechanism  formation  to  the  point  of  collapse  (2).  The  maximum  load  predicted 
by  BEAM  will  therefore  approximate  the  collapse  load  of  the  frame. 


0  20  40  60  80  100  120 

Horizontal  displacement,  mm 

Figure  10.  An  enlurgcnu'nt  of  the  initial  portion  of  the  load-displacement  curve  for 
Frame  4,  comparing  elastic,  rigid-plastic  and  Step-B}/-Step  anali/ses  with  experimental 
results.  The  combined  model  of  the  load-displacement  curve  is  shown  by  the  bold  line.  The 
rigid-plastic  model  assumes  that  the  plastic  hinges  form  at  the  corners  of  the  frame  and  its 
fixed  supports. 


A  compatibility  check  was  also  performed  to  check  the  reliability  of  the  joint 
and  member  rotations  predicted  by  BEAM.  These  rotations  were  obtained  from 
the  output  file  in  Section  3  of  the  Appendix,  converted  into  degrees,  and 
compared  with  the  axis  rotations  which  were  the  rotations  of  the  element 
centrelines  from  their  original  positions.  Fig.  11.  These  comparisons  were  made  at 
the  point  of  mechanism  formation.  The  following  compatibility  check  shows  that 
all  these  rotations  are  consistent,  where  AR,  JR,  and  MR,  refer  to  axis,  joint  and 
member  rotations,  respectively: 

element  1,  node  2  AR  =JR-rMR  -2.16°  -  -1.17°  -r  (-1.17°) 

element  2,  node  2  AR=MR-JR  0°  =  -1.17°  -  (-1.17°) 


element  2,  ntxle  3  AR=MR-JR  0°  =  -1.16°  -  (-1.16°) 
element  3.  ncxle  3  AR=JR+MR  -2.16°  ~  -1.16°  +  (-1.16°) 


♦; 


•  • 


_  • 
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38.83 


38.74 


Figure  11:  Axis  rotations  at  the  point  of  mechanism  formation.  Dimensions  in  mm. 


The  rigid-plastic  model  that  was  depicted  in  Fig.  9,  predicted  a  plastic  hinge 
rotation  of  2.34‘^  for  a  39  mm  frame  deflection  at  a  load  of  12.48  kN.  It  is 
interesting  to  note  that  this  rotation  is  equivalent  to  the  axis  rotation  obtained  by 
adding  the  joint  and  member  rotations  at  element  1,  node  2.  The  rigid-plastic 
model  predicted  collapse  to  occur  at  a  plastic  hinge  rotation  of  5  95“,  which 
corresponds  to  a  98  mm  deflection  of  the  frame  at  a  load  of  12  33  k.\ 

In  general,  there  is  good  correspondence  between  the  predictions  of  the 
combined  model  and  the  experimental  results.  While  the  stiffness  of  the  frame  is 
overpredicted  prior  to  collapse,  this  is  because  its  joints  and,  especially  the  fixity 
at  Its  base,  are  considered  to  be  perfectly  rigid.  Furthermore  the  elastic  and  Step- 
By-Step  analyses  are  only  first-order  models  and  do  not  therefore  account  for  the 
changing  frame  geometry.  The  differences  such  as  there  are  between  the  predicted 
and  experimental  load -displacement  curves  beyond  collapse  are  mainly  due  to  the 
design  of  the  rigid  )oints  of  the  frame  which  caused  the  plastic  hinges  to  develop 
in  slightly  different  locations  to  those  predicted  by  BEA.M,  i.e.  at  the  edges  of  the 
rigid  joints.  This  was  to  prevent  joint  failure  such  as  had  occurred  in  earlier  frame 
tests  (1),(2). 

Overall  it  can  be  seen  that  BEAM  provides  a  valuable  adjunct  to  separate  rigid- 
plastic  modelling.  As  well  as  predicting  the  locations  of  the  plastic  hinges,  BEAM 
also  provides  an  accurate  estimate  of  the  collapse  load  of  the  frame,  i.e.  12.18  kN, 
compared  to  12,20  kN  found  experimentally.  This  is  because  BEAM  should 
usually  satisfy  the  mechanism,  yield  and  equilibrium  conditions.  Moreover,  if  the 
load-deflection  curve  produced  by  BEAM  is  combined  with  a  rigid-plastic  model 
of  the  collapsing  frame  then  the  definition  of  the  elasto-plastic  transition  prior  to 
collapse  is  improved  compared  to  simpler  modelling  methods  (1),(2).  Figures  9 
and  10  show  that  simply  sketching  a  curve  by  hand  from  the  point  of  mechanism 
formation  (predicted  by  BEAM)  to  the  actual  point  of  collapse  (predicted  by  the 
separate  rigid-plastic  analysis)  completes  the  prediction  of  the  frame’s  load- 
displacement  curve  (1),(2). 
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5.  Conclusion 


Program  BEAM  will  determine  the  deformation  of  a  thin-walled  frame  up  to  the  point  at 
which  sufficient  plastic  hinges  form  to  transform  the  frame  into  a  mechanism,  i.e.  the 
point  of  mechanism  formation.  This  allows  the  locations  of  the  plastic  hinges  ct>mposing 
the  collapse  mechanism  of  the  frame  to  be  predicted  along  with  an  estimate  ot  its 
maximum  load.  When  BEAM  is  used  in  conjunction  with  a  large-deflection  rigid-plastic 
mcxlelling  technique,  the  U>ad-deflection  cur\’e  of  a  thin-walled  frame  can  be  determined 
beyond  the  jxiint  of  mechanism  formation  to  where  drcx>ping  livad-displacement 
behaviour  is  exhibited  at  large  deflections 


( 


♦ 


•  { 
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Appendix  1 

Program  BEAM:  Input  File  Layout 


PROGRAM  BEAM 


Th«  following  t«xt  d««crib«a  th«  input  data  layout: 
CARD  HO. -daacription  COMMAMO  SYNTAX 


1-titla 
2*- control 

3 - coordinataa 

4- iiatarials 

5- gaonatric  prop. 

6- connactivity 

7- iBaabar  act  Iona 

8- 

9*-nodal  actions 

lo¬ 
ll -aattlaaant  a 
12- 

13  -diaplacaitanta 


9n  COMMAMO  SYNTAX 

TITLE 

NNODE « NNODSL , NMAT , NGBO , MELEM , NBD 

I.X. Y.Z 

I.E.PR.G 

Z<AA. J,I2,I3«MP 

Z  «  NX , N J , NX , MT YP  <  STYP , BTYP 

M 

N,  PEAl,  PBA2 . PEA12 

M 

N, P1,P2 . P12 

N 

O ,  VALl , VAL2 , VAL3 . VAL4 , VAL5 , VAL6 
N.  DOPl .  DOP2  ,  OOP 3  .  DOP4  .  DOP5,  DOP6 


Whara:  TITLE 

NNODE 

NMODBL 

NMAT 

NGBO 

NELEM 

NBD 


probXan  titia  (60  characters  loog> 
nuoibar  of  structural  nodes 

total  nuAbar  of  nodes  listed  (incl.  k  nodes) 
nuisbar  of  aatarial  types 
nuAbar  of  gaoAatric  sactions 
nuAber  of  alaAenta 

ouAbar  of  nodal  diaplacenent  conditions 

Node  No.  ( -va  for  k  nodes) 

X  coordinate  of  node  <a.g.  aatras> 

Y  coordinate  of  node  <a.g.  natrea> 

Z  coordinate  of  node  <a.g.  Aatras> 

M.B.  The  nodes  should  be  in  sequential  order 


c 

I 

- 

Material  No. 

c 

E 

(MAT (I, 1)1 

Young's  modulus  <a.g.  kNm~^> 

c 

PR 

IMAT(I,2)] 

Poisson's  ratio 

c 

c 

c 

G 

[MAT(I,3)] 

Shear  Modulus  <a.g.  kHa~  >  [Optional!  If 
to  xaro,  G  calculated  from  E/2<1«^PR)] 

c 

c 

I 

Geometric  Property  No. 

c 

AA 

[GEO (I, 1)1 

Axial  Area  <a.g.  a^> 

c 

J 

[OEO(I,2)l 

Moment  of  Inertia  abt.  1  axis  <a.g.  m^> 

c 

12 

(GEO (I, 3)1 

**  2  axis  «a.g.  m^>^ 

c 

13 

(OEO(I,4)J 

*  3  axis  <e.g.  m^> 

c 

c 

MP 

[GEO(I,5)] 

Failure  Moment  <e.g.  kNm» 

c 

c 

I 

Eleaent  No . 

c 

MI 

C ICON (1,1)1  - 

Global  node  number  I 

c 

NJ 

[ ICON (1,2)1  - 

-  J 

c 

NX 

[ ICON (1,3)1  ' 

**  K  ( ♦ve  in  ♦ve  quad. 

c 

MTYP  ( ICON (1,4)1  - 

Material  No.  -ve  in  -ve  quad) 

c 

STYP  (ICON (I, 5)1  - 

Geometric  Property  No. 

c 

c 

BTYP  (ICON (I, 6)1  - 

Eleaent  Type  (Set  to  zero  as  not  used] 

c 

c 

c 

M 

-  no . 

of  elements  subjected  to  Member  Actions 

<e.g.  kN  or  kNa> 

-  alaaent  no. 


It  Dor  1 


e 

c 

e 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C  OMLY  FOR 
C  ELASTIC 
C  ANALYSIS 
C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C*********' 


riAi  (A(1)] 

I 

I 

rSAlE  CA(12)] 


-  FijMd  Bad  Aetioa  at  ali 
(local  coord)  I 

I 

•  rlxad  Bad  Aetioa  at  alcAcat  OOP  13 
(local  coord) 


M 

•  no.  of  nodaa  aubjactad  to  Nodal  Actiona 

<a.g.  ktt  or  )Kna> 

¥ 

-  noda  no. 

PI 

CA(1)] 

1 

-  Nodal  Action  at  noda  OOP  1 

1 

P6 

1 

tA<C)] 

1 

-  Nodal  Action  at  noda  OOP  € 

¥  -  no.  node*  cubjaetod  to  aattlaaMnta/ 

apaeifiad  dlaplaeaaanta 

N,VAL1«VAL3.VAL3,VAL4.VALS.VAL6  noda  no..  (IDOP1-IDOP6] 

If  VALna'val*.  thia  OOP  ia  aat  to  'val*. 
tha  aattlaaant  .<a.g.  aatraa> 

If  VALna‘0‘«  no  valua  ia  aat  for  thia  OOP 
ag.  y  aattlaaant  of  0.05a  :  0.0.0.05.0.0.0,0 


N.0OP1.OOP3.0OP3.0OP4. OOPS. OOPS  noda  no..  (XDOP1-IDOP6] 
If  DOPn«‘0‘.  thia  OOP  ia  aat  to 
saro  dlaplacaaant 

Zf  OOPna'l*.  tha  aovaaant  of  thia  OOP 
la  poaaibla  and  will  ba  datarainad 
ag.  30  flxad  and  t  M. 0. 0. 0, 0 . 0 . 0 
ag.  30  aiapla  aupport  :  M. 0. 0. 0. 1. 1.  1 
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Appendix  2 

Program  BEAM:  Input  File 


2D  PRAMB 
4,7, 1,1, 3,2 
1,0, 0,0 
2,0,1.025,0 
3,1.180,1.025,0 
4,1.180,0,0 
5, -0.2, 0.5,0 
6,0  5, 1.5,0 
7, 1.5, 0.5,0 
1,20686,0.3,0 

1, 3848-6,0.22128-3, 1477128-12, 1477128-12, 3.12 
1, 1,2,5, 1,1,0 
2, 2, 3, 6, 1.1,0 
3. 3, 4, 7. 1,1,0 
0 
1 

2,5, 0, 0, 0,0,0 
0 

1,  0,  0,  0,  0,0,0 

4, 0, 0, 0,  0,0,0 


•  • 
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Appendix  3 

Program  BEAM:  Output  File 


2D  nUME 


CONTROti  DATA 


mMBBR  OF  MOORS  ZN  STRUCTURE  a  4 
TOTAL  MUMBIR  OF  MOORS  (ZMCL.  K  MOORS)  «  7 
NUMBER  OF  MATRRIA1«S  -  1 
NUMBER  OF  OROMRTRZC  SRCTZONS  a  1 
NUMBER  OF  BRAN  ELEMENTS  a  3 
NUMBER  OF  BOUNDARY  CONDITIONS  a  2 


NODAL  COOROZMATES  TABLE 


NODE 

X 

Y 

z 

1 

0.0000 

0.0000 

0.0000 

2 

0.0000 

1.0250 

0.0000 

3 

i.iaoo 

1.0250 

0.0000 

4 

1.1800 

0.0000 

0.0000 

5 

-0.20 

0.5000 

0.0000 

6 

o.sooo 

1.5000 

0.0000 

7 

l.SOOO 

0.5000 

0.0000 

MATERIAL 

PROPERTIES 

MATERIAL 

YOUNG  S 

POISSON  S 

SHEAR 

MODULUS 

RATIO 

MODULUS 

1 

0.2060E«-09 

0 

.  3000 

0 

.7923E*0e 

GEOMETRIC 

PROPERTIES 

SECTION 

AX.  AR.  J 

12 

13 

MF 

1  0.3a40E-03  0.2212E-03 

0. 1477E 

-06 

0.1477E-06  0 

.3120E4-01 

ELEMENT  CONNECTIVITY 

TABLE 

ELEMENT 

NODEl  NODE2 

K  1 

MATERIAL  SECTION 

ELEMENT  TYPE 

1 

1 

2 

5 

1  1 

0 

2 

2 

3 

6 

1  1 

0 

3 

3 

4 

7 

1  1 

0 

APPLIED  NODAL  ACTIONS 

IN  TERMS  OF  GLOBAL  AXES 

JOINT  NO. 

X  FORCE  Y 

FORCE 

Z  FORCE 

X 

MOMENT  Y  MOMENT  Z  MOMENT 

2 

5.00 

0.00 

0.00 

o 

o 

o 

o 

o 

o 

o 

o 

o 
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KZMtMXTZC  BOUNDARY  CONDZTZONS 


NODS  DOPX  DOP2  DOP3  DOP4  OOPS  OOP6 

1  0.  0.  0.  0.  0.  0.  ^ 

4  0.  0.  0.  0.  0.  0. 

•••••••*•  QP  ihput  data  •••••••*• 

PXJ^TZC  HZNOE  NO.  I  IS  PORMEO  IN  ELEMENT  I  % 

AT  NODE  1  WHEN  LOAD  PACTOR>  2.09310 

CUMULATIVE  DISPLACEMENTS 


TRANSLATIONS  ROTATIONS 


JOINT 

X 

Y 

Z 

X 

Y 

z 

1 

O.OOOOB^OO 

O.OOOOE«-O0 

O.OOOOE«0O 

O.OOOOB^OO 

O.0000E«^00 

O.OOOOB4^00 

3 

0.3396B-01 

0.4937B‘04 

O.OOOOE^OO 

O.OOOOE*00 

0  .  OOOOEt^OO 

-0. 1465E-01 

3 

0.3388E-01 

•0.4937E-04 

O.OOOOEfOO 

O.OOOOE«-O0 

O.OOOOE4>00 

-0.1457E-01 

4 

O.OOOOE^OO 

O.OOOOE^OO 

0.0000E4-00 

O.OOOOE^OO 

O.OOOOBfOO 

O.OOOOB^OO 

CUMULATIVE  STRESS*RESULTANTS 


•  • 


HEM  MODE  AXIAL  X  MOMENT  Y  MOMENT  Z  MOMENT  PAZLURE  MOMENT 


1  1  •0.3810e<»01  0.0000E«-00  O.OOOOE4>00  0.3120E4>01  0.3120B4>01 

1  2  0.3810E4'01  O.OOOOE>00  O.OOOOBfOO  0.3250Xf01  0.3120B^01 

3  2  0.5336E4-01  O.OOOOE^-OO  0.0000E«>O0  -0.3350Bf01  0.3130B«-01 

3  3  -0.5326E4-01  O.OOOOE>00  0.0000E4-00  -0.3346E4>01  0.3130E«^01 

3  3  0.3810E>01  O.OOOOE^OO  O.OOOOE  +  00  0.3346E«>01  0.3120B4-01 

3  4  >0.3810E^01  O.OOOOE^OO  0.0000E«00  0.3111E«01  0.3130B^01 

CUMULATIVE  MEMBER  DISPLACEMENTS 


MEM  NODE  AXIAL  X  ROT  Y  ROT  Z  ROT 


1  1  O.OOOOB+OO  O.OOOOB^-OO  O.OOOOE^OO  O.OOOOB^OO 

1  3  0.4937B'04  O.OOOOE+OO  0.0000E«00  -0.146SB-01 

3  3  0.3396B-01  O.OOOOB4-00  O.OOOOE^OO  -0.1465E-01 

3  3  0.3388E-01  O.OOOOB4-00  0.0000E«-00  -0.1457E-01 

3  3  0.4937B-04  0.0000E*00  O.OOOOEtOO  -0.1457E-01 

3  4  O.OOOOEtOO  0.0000B«-00  O.OOOOB^OO  O.OOOOB-^00 


3 


rLMTZC  aZNOB  MO.  2  Z8  rOMKD  ZN  BLBHBMT 
AT  MODI  4  WHBN  LOAD  PACTOA«  2.09682 


CUMULATZVB  DZ8PLACBNBNT8 


TlUUf8IATZONS  ROTATZOM8 


JOINT 

X 

Y 

Z 

X 

Y 

X 

1 

O.OOOOB^OO 

0.00008*^00 

0.00008*00 

0.00008*00 

0.00008*00 

0.00008*00 

2 

0.2303B-01 

0.49488-04 

0.00008*00 

O.OOOOE*00 

0.00008*00 

-0.14678-01 

3 

0.22958-01 

-0.4948B-04 

0.0000E*00 

0.00008*00 

0.00008*00 

-0.14628-01 

4 

0.0000B«^00 

0.00008*00 

0.0000B*00 

0.00008*00 

0.00008*00 

0.00008*00 

cumu««a;x’/b  strj^sS'Resultamts 


MEM  NODE  AXZAL  X  MOMENT  ¥  MOMENT  Z  MOMENT  PAZLORB  MOMENT 


I  -I  -0.3819E*01  0.0000£«00  O.OOOOE^OO  0.3120B^01  0.3120B«>01 

1  2  0.38I9E1-01  O.OOOOBtOO  0.0000E*00  0.2254B+01  0.3120E4-01 

2  2  0.5241E^01  O.OOOOB^OO  0.0000B«00  *0.2254B«>01  0.3120X«01 

2  3  -0.S241E  +  01  O.OOOOE^-OO  0.00O0E«0O  -0.22S2B*01  0.3120E4-01 

3  3  0.3ei9Et01  0.0000E*00  O.OOOOE  +  00  0.2252E«‘01  0.3120E«-01 

3  4  -0.38X9E»01  O.OOOOB^-OO  O.OOOOE+00  0.31208i-01  0.3120B4-01 

CUMULATZVB  MEMBER  OZSPLACBNENTS 


MEM  NODE  AXZAL  X  ROT  Y  ROT  Z  ROT 


1  -1  O.OOOOBoOO  O.OOOOEfOO  O.OOOOE^OO  0.0000£«>00 

1  2  0.4948E-04  O.0OOOE4>OO  O.OOOOE^-OO  •0.I467E-01 

2  2  0.3303E-01.  O.OOOOE^OO  O.OOOOE^OO  •0.14678-01 

2  3  W.2295S-01  O.COCOCi>00  0.00O0Ei>00  -0.1462Z  '' ' 

3  3  0.4948E'04  O.OOOOEt-OO  O.OOOOE^OO  -0.1462E-0X 

3  4  O.OOOOE^OO  O.OOOOE^OO  O.OOOOE^-OO  O.OOOOE^OO 


PLASTIC  HINGE  NO.  3  IS  FORMED  IN  ELEMENT  2 
AT  NODE  2  WHEN  LOAD  PACTORs  2.43456 


CUMULATIVE  DISPLACEMENTS 


TRANSLATIONS 

ROTATIONS 

JOINT  X 

Y 

Z 

X 

Y 

z 

1  O.0000E*O0 

0.0000E*00 

0.OOOOE*OO 

0.0000E*00 

0.00008+00 

0.00008*00 

2  0.38768-01 

0.6849E-04 

O.OOOOE*00 

O.OOOOE*00 

O.OOOOE+OO 

-0.20308-01 

3  0.3867E-01 

-0.6849E-04 

0.00008*00 

0.00008*00 

O.OOOOR+00 

-0.20248-01 

4  0.0000E*00 

0.00008*00 

0.00008*00 

o.ooooe*oo 

0.00008*00 

0.00008*00 

CUMOLATIVE  STRESS ’RESULTANTS 


MEM 

NODE 

AXIAL 

X  MOMENT 

Y  MOMENT 

Z  MC»(BMT 

FAILURE  MOMENT 

1 

-1 

-0.52a6E4-01 

O.OOOOEt-00 

O.OOOOE+00 

0.3120E«‘01 

0.3120B4>01 

1 

2 

0. 5266B4'01 

0.O000B-»^0O 

0 . 0000E«00 

0.3120Er01 

0. 3120E4-01 

2 

2 

0.6085E«01 

O.OOOOEf^OO 

O.OOOOB^OO 

-0.3120E*01 

0.3130E4‘01 

3 

3 

-0.6085E^0I 

O.OOOOE^OO 

O.O0OOE»OO 

-0.3117B+01 

0. 3120E4-01 

3 

3 

0.5286E«-01 

O.OOOOBi-00 

O.OOOOC+00 

0.3117B»01 

0.3130E4-01 

3 

-4 

-0.5286E>01 

O.OOOOE+00 

O.OOOOEi-OO 

0.3120E«^01 

0.3120E4-01 

CUMULATIVE  MEMBER  DISPLACEMENTS 

MEM 

NODE 

AXIAL 

X  ROT 

Y  ROT 

Z  ROT 

1 

-1 

O.OOOOE^OO 

0.0000E4-00 

0.O000E*0O 

O.OOOOB^OO 

1 

2 

0.6849E-04 

O.OOOOB4-00 

0 . 00O0E*O0 

-0.2030E-01 

2 

2 

0.3876E-01 

O.OOOOE^OO 

O.OOOOE«00 

-Q.2030E-01 

2 

3 

0.3667E-01 

0.0000E4-00 

O.OOOOEf^OO 

’0.2024C-D1 

3 

3 

0.6849B-04 

O.OOOOBt-OO 

O.OCOOE^OO 

-0.2024E-01 

3 

-4 

O.OOOOE^OO 

O.OOOOE^OO 

O.OOOOE^OO 

O.OOOOE.OO 

PLASTIC  HINOE  NO.  4  IS  FORMED  IN  ELEMENT  2 
AT  NODE  3  WHEN  LOAD  FACTOR*  2.43513 


COKULATIVE  DISPLACEMENTS 


TRANSLATIONS  ROTATIONS 


JOINT 

X 

Y 

Z 

X 

Y 

z 

1 

O.OOOOE«-O0 

O.OOOOE^OO 

O.OOOOE^OO 

O.OOOOE^OO 

O.OOOOB-t-OO 

0.0000E4-00 

3 

0. 3883E-01 

0.6852E-04 

O.OOOOE^OO 

0.OOOOE»OO 

O.OOOOE4-00 

-0.2037E*01 

3 

0. 3874E-01 

-0.6852E-04 

O.OOOOB^OO 

0.0000B4-OO 

O.OOOOB4-00 

’0.3038B-01 

4 

O.OOOOE^OO 

O.OOOOE^OO 

0.0000E«-00 

O.OOOOB^OO 

O.OOOOB^OO 

0.0000B-*>00 

CUMULATIVE  STRESS’RESULTANTS 


MEM 

NODE 

AXIAL 

X  MOMENT 

Y  MOMENT 

Z  MOMENT 

FAILURE  MOMENT 

1 

-1 

-0.5288E4-0X 

O.OOOOE+00 

O.OOOOE+OO 

0.3130E+01 

0.3130E+01 

1 

2 

0.5388E1-01 

O.OOOOE+00 

O.OOOOE+OO 

0.3120E+01 

0.3130B+01 

3 

-2 

0.60B8E4-01 

O.OOOOB+00 

O.OOOOE+00 

’0.3130E+01 

0.3120E+01 

3 

3 

-0.6088Bt01 

O.OOOOB+OO 

O.OOOOB+OO 

-0.3130E+01 

0.3120E+01 

3 

3 

0.53888+01 

O.OOOOE+OO 

O.OOOOB+00 

0.3130E+01 

0.3120B+01 

3 

>4 

’0.5288B+01 

O.OOOOB+OO 

O.OOOOB+00 

0.3I30B+01 

0.3130B+01 

COMUIATXVS  OXatLACIMBIITS 


MM  MODS  AXXAL  X  XOT  Y  XOT  S  ROT 


1  -1  O.OOOOB^OO  O.OOOOB^OO  0.0000s«>00  O.OOOOB^OO 

1  a  0.68S3B-04  0.0000B«00  O.OOOOI^OO  -0.30371-01 

3  -3  0.3t83B-01  0.0000B«^00  O.OOOOB^OO  -0.3037B-01 

3  3  0.3874B-01  O.OOOOl^OO  O.OOOOl^OO  -0.3038X-01 


3  3  0.68S3B-04  O.OOOOB^OO  O.OOOOC4>0O  -0.3038B-01 
3  -4  O.OOOOB^OO  O.OOOOK^OO  O.OOOOB^OO  O.OOOOX^OO 


DBPOXMATXONS  LAXOBX  THAN  0.110B^03  IN  CYCLE  MO.  5 


CUNULATIVB  DISPLACBMBMTS 


# 

TllANSLATXOHS  ROTAYIOHS 

JOINT  X  Y  Z  X  Y  Z 

1  O.OOOOE^OO  0.0000X4-00  0.0000X400  0.0000X400  0.0000X400  0.0000X400 

2  0.3883X-01  0.66528-04  O.OOOOX4OO  O.OOOOB4OO  O.OOOOB4OO  -0.a037X-01 

3  0.3874E-01  -0.6852X-04  O.OOOOX4OO  O.OOOOE4OO  O.OOOOB4OO  -0.2028X-01 

4  0.0000X400  0.0000X400  O.OOOOX4OO  O.OOOOX4OO  O.OOOOX4OO  0. 0000X400  0 

CUMULATIVE  STRESS-RESULTANTS 


MEM  NODE  AXIAL  X  MOMENT  Y  MOMENT  Z  MOMENT  FAILURE  MOMENT 


1  -1  -0.S28eE401  O.OOOOE4OO  O.OOOOE4OO  O.312OE4OI  0.3120X401 

1  2  O.S288E401  O.OOOOE4OO  O.OOOOB4OO  0.3120B401  0.3120B401 


•  • 


2  -2  0.6088E40i  O.OOOOE4OO  O.OOOOE4OO  -0.3120X401  0.3120X401 

2  -3  -O.6O88E4OI  O.OOOOE4OO  O.OOOOB4OO  -0.3120X401  0.3120X401 


3  3  0.52868401  O.OOOOB4OO  O.OOOOE4OO  0.3120X401  0.3120X401 

3  -4  -O.S28eE401  O.OOOOE4OO  O.OOOOE4OO  0.3120E401  0.3130B401 

CUMULATIVE  MEMBER  DISPLACEMENTS 


MEM  NODE  AXIAL  X  ROT  Y  ROT  Z  ROT 


1  -1  O.OOOOE4OO  O.OOOOC4OO  O.OOOOE4OO  0.0000E400 

1  2  0.6852E-04  O.OOOOE4OO  O.OOOOE4OO  -0.3037E-01 

2  -2  0.36838-01  O.OOOOE4OO  O.OOOOE4OO  -0.2037E-0I 

2  -3  0.3874E-01  O.OOOOE4OO  O.OOOOE4OO  -0.2028B-01 

3  3  0.6852E-04  O.OOOOE4OO  O.OOOOE4OO  -0.3028C-01 

3  -4  0.00006400  0.00006400  O.OOOOE4OO  O.OOOOS4OO 


UNPACTORBD  DISPLACEMENTS  AFTER  LAST  SOLUTION 

• 

TRANSLATIONS 

ROTATIONS 

JOINT  X 

Y 

Z 

X 

Y 

z 

1  0.00006400 

2  -0.4993E412 

3  -0.4993E412 

4  O.OOOOE4OO 

0.00006400 

-0.27886-22 

-0.54736-23 

O.OOOO64OO 

0.0000X400 

0.00006400 

O.OOOOE4OO 

0.00006400 

O.OOOOE4OO 

O.OOOO64OO 

O.OOOOB4OO 

O.OOOOX4OO 

O.OOOOE4OO 

0.0000X400 

O.0000E400 

0.0000X400 

0.0000X400 

0.4871X412 

0.4871X412 

O.OOOOE4OO 

• 
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